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Rewetting global wetlands effectively reduces 
major greenhouse gas emissions 


Junyu Zou®', Alan D. Ziegler ®?, Deliang Chen®?, Gavin McNicol ©‘, Philippe Ciais*°, 
Xin Jiang’, Chunmiao Zheng ®"”, Jie Wu", Jin Wu?"®, Ziyu Lin"’, Xinyue He®"", Lee E. Brown ®”, 
Joseph Holden®”, Zuotai Zhang’, Sorain J. Ramchunder™, Anping Chen®™ and Zhenzhong Zeng®'™ 


Carbon and nitrogen losses from degraded wetlands and methane emissions from flooded wetlands are both important sources 
of greenhouse gas emissions. However, the net-exchange dependence on hydrothermal conditions and wetland integrity 
remains unclear. Using a global-scale in situ database on net greenhouse gas exchanges, we show diverse hydrology-influenced 
emission patterns in CO,, CH, and N,O. We find that total CO,-equivalent emissions from wetlands are kept to a minimum when 
the water table is near the surface. By contrast, greenhouse gas exchange rates peak in flooded and drained conditions. By 
extrapolating the current trajectory of degradation, we estimate that between 2021 and 2100, wetlands could result in green- 
house gas emissions equivalent to around 408 gigatons of CO,. However, rewetting wetlands could reduce these emissions 
such that the radiative forcing caused by CH, and N,O is fully compensated by CO, uptake. As wetland greenhouse gas bud- 
gets are highly sensitive to changes in wetland area, the resulting impact on climate from wetlands will depend on the balance 


between future degradation and restoration. 


etlands have continuously accumulated organic car- 

bon since the Last Glacial Maximum’, forming a dense 

carbon pool that stores over one-third of global soil 
organic carbon in only 6% of the total land area’. Alarmingly, since 
the Industrial Revolution, more than half of wetlands have been 
degraded by anthropogenic activities, including drainage, defores- 
tation, afforestation, agricultural expansion, urbanization and cli- 
mate change**. A phenomenon commonly associated with wetland 
degradation is the lowering of the water table, which exposes car- 
bon pools above the water table to decomposition and releases CO, 
while simultaneously altering the natural exchange of other green- 
house gases (GHGs) including methane (CH,) and nitrous oxide 
(N,O)°-". Given that degraded wetlands are important sources of 
GHG emissions to the atmosphere''"’, there is a critical need not 
only to determine the impact of widespread wetland degradation on 
GHG exchanges, but also to assess the potential for wetland restora- 
tion in reducing GHG emissions. 

Fluxes of CO,, CH, and N,O in wetlands are mediated by the 
water level relative to the surface”’>’*. Electron acceptor limitation 
in the soils of saturated wetlands is favourable to the production of 
methane'”'’, forming the world’s largest natural source of CH, emis- 
sions'*. Meanwhile, substantial CO, is consumed through photosyn- 
thesis by wetland vegetation; and the anaerobic conditions reduce 
the decomposition of organic carbon, generating a major terrestrial 
carbon sink over long timescales'*'°°. By contrast, the degradation 
of wetlands, either through drainage or desiccation, exposes stored 
organic matter to aerobic decomposition, resulting in the emission 


of CO, and N,O to the atmosphere'"'*"", Previously, the relation- 
ship between GHG emissions and the water table in wetlands has 
been examined for only one or two GHGs”, and for one single wet- 
land category (for example, managed peatland’***) or various wet- 
land categories at local scales''’>’’. Lacking is a global assessment 
involving the three main GHGs to guide the development of effec- 
tive climate change mitigation strategies and to inform the potential 
to restore the functioning of wetland ecosystems across moisture 
and temperature regimes worldwide. 

We address this issue by building a global database containing 
in situ observations of exchange rates of GHGs for wetlands, drawn 
from 3,704 site-year records (Extended Data Fig. 1; refer to Methods 
for details). Each record contains data on wetland environmental 
conditions and flux information that allows quantitative assessments 
of the net fluxes of CO,, CH, and N,O under various wetness condi- 
tions. The reported details of the environmental conditions from each 
site-year record allow us to perform a new multi-gas assessment for 
a variety of wetland types and moisture regimes worldwide. In doing 
so, we classify the wetness condition of each site-year in the growing 
season into one of six categories that are related to the water table/ 
level (WTL) depth below (negative number) and above (positive) 
the surface: WTL-3, <-70cm; WTL-2, -70 cm < WTL-2 <-50cm; 
WTL-1, -50cm<WTL-1<-30cm; WTLO, -30cm<WTLO<-5 
cm; WTLI1, -5cm<WTL1<40cm; and WTL2, >40cm. We 
illustrate differences related to temperature regimes by assessing 
responses across three climate zones (boreal, temperate and tropical; 
defined by thresholds of multiyear-average surface air temperature). 
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Fig. 1] WTL effects on global wetland NEE and total GHG emissions. a, Relationship between NEE and sum of three GHG (CO,, CH, and N,O) net fluxes 
in different WTLs, drawn from 174 site-year records that reported three GHGs. b, NEE, CH,, N,O and sum of three net fluxes for different WTL conditions. 
c-f, Total (¢) and individual (CO, (d), CH, (e), NO (f)) GHG fluxes for the six different WTLs considered. Points in each box are sampled from the original 
dataset (3,672 site-year records total) with 1,000 bootstraps. Different letters in the boxes indicate significant differences (P< 0.01) between various 
WTLs based on nonparametric Wilcoxon signed-rank tests. Bold vertical lines show the median; boxes indicate the middle two quartiles; horizontal lines 
indicate the non-outlier range. Note that x axes have been truncated for enhanced readability. 


The nonlinear hydrothermal influence on GHG exchange 
We call the CO, net-exchange flux at the water/land-atmosphere 
interface as net ecosystem exchange (NEE) (equation (1); positive/ 
negative values indicate GHG sources/sinks). Through establishing 
relationships between NEE and the total GHG flux (sum of CO,, 
CH, and N,O in CO,-equivalent (CO,eq)) based on the records 
containing complete data, we observe the following: (1) almost all 
(173 of 174) records show the total GHG flux values exceed NEE 
for the same site-year, and (2) the differences between the NEE 
and total GHG flux are highly dependent on moisture conditions 
(Fig. la,b). These results indicate clearly that the wetness regula- 
tion pattern for CH, or N,O emissions is different from that of CO, 
(ref. *”). By mapping the distribution of wetness control for the three 
GHGs for various types of wetlands (inter alia bogs, fens, marshes, 
swamps, floodplain and water bodies; Methods), we identify non- 
linear (parabolic) exchange patterns for NEE and the sum of GHGs 
and opposing monotonic patterns for CH, and N,O (Fig. 1c-f and 
Extended Data Figs. 2a,b, 3b and 4a,b). Maximum emission of 
CH, occurs when flooded wetlands have water levels well above 
the soil surface (WTL2; water level> 40cm), and the minimum 
occurs when the water table is well below the surface of the wetland 
(WTL-2 and WTL-3; water table <—50 cm), indicative of a drained 
or desiccated state. By contrast, the highest emissions of N,O occur 
during dry conditions (WTL-3), and the lowest occur in flooded 
conditions (WTL2). Emissions of CO, exhibit relative extremes for 
both high water-level and low water-table conditions (Fig. 1d and 
Extended Data Fig. 4a). 

These hydrology-dependent emission patterns are in line with the 
expectation that CH, is produced in the anaerobic conditions that are 
associated with waterlogged soils*'”*. Lower emissions of N,O occur 


during flooded conditions because facultative anaerobic denitrifying 
bacteria reduce N,O to N, in the oxygen-depleted water column’. 
As the water table falls below the soil surface, aerobic decomposi- 
tion of organic matter results in an increase in CO, emissions'!">?!. 
As expected, the highest CO, emission was observed in WTL-3, 
where the water table is lower than 70cm (Fig. 1d and Extended Data 
Fig. 4a). The relatively high CO, emissions observed in wetlands 
under flood conditions (WTL2) are probably driven by the lateral 
movement of organic matter across the landscape, the leaching of 
organic carbon into a dissolved state and subsequent oxidation by 
heterotrophs**’. Finally, the lowest emission of the sum of all three 
GHGs (CO,, CH, and N,O) occurs when the water table is near the 
ground surface (WTLO; ranging from -30cm to -5cm), with near 
zero emissions. This nonlinear wetness pattern is in agreement with 
the recent report by ref. *, which studied only selected peatland sites 
and did not incorporate N,O (Extended Data Fig. 2a,b). 

The parabolic pattern in GHG flux varies across boreal, tem- 
perate and tropical regimes (Fig. 2 and Extended Data Fig. 3b), 
demonstrating that climate influences the dependency of wet- 
land GHG emissions on the wetness regime. The net GHG flux in 
each temperature regime tends to be approximately neutral when 
the water table is near the ground surface, although it requires 
higher WTL (-5 to 40cm) to approach GHG equilibrium for 
the tropical sites (Fig. 2c and Extended Data Table 1). Using the 
empirical GHG exchange rates in the WTL2 group for different 
temperature regimes (Methods), we estimate the annual GHG 
emissions from global water bodies (lakes and reservoirs) to be 
1.0 GtCO, yr! and 127.5 TgCH, yr", which is similar to the previous 
reports**’ of 1.2GtCO,yr! and 175.2TgCH, yr (Supplementary 
Fig. 1). In addition, we estimate the CH, emissions from natural 
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Fig. 2 | Nonlinear hydrothermal influence on GHG exchange. a-c, Dependency of GHG emissions in boreal (long-term average air temperature <4 °C) 


(a), temperate (long-term average air temperature 4-17 °C) (b) and tropical 


(long-term average air temperature >17 °C) (¢) regions on WTL and 


climate. d, The ‘mean’ groups are calculated from equilateral weighted averages in each climate regime. Dots and shadows represent mean +1.96 SEs. 
e-h, Contribution ratios of NEE, CH, and N,O to the sum of three GHG net fluxes in the three climatic regions (boreal (e), temperate (f) and tropical (g)) 


and the mean Ch). 


freshwater wetlands to be 144.4Tg yr", corresponding well with the 
148.6 TgCH, reported by ref. '*. Furthermore, the GHG emissions 
per area for wetlands with a low water table (WTL <-70 cm) are 19.7 
and 11.2tCO,eqha"' yr" for boreal and temperate regimes, respec- 
tively. These values are consistent with both ref. * (17.60 tCO,eq 
ha"'yr7 in the case of -70cm water level) and ref. '’, who deter- 
mined the drainage-related GHG emission rates for boreal and 
temperate zones (16.1tCO,eqha™yr-'). These validated empirical 
values in the nonlinear relationship between WTL and GHG emis- 
sions then represent a novel opportunity to assess the global GHG 
emissions resulting from wetland degradation. 


GHG emissions from wetland degradation 

We evaluate past and future scenarios of wetland degradation by 
integrating the natural Wetland Extent Trends (WET) index*? 
(Supplementary Fig. 2) with the Global Lakes and Wetlands 
Database** (GLWD). We assess the historical emissions from 
degraded wetlands at the global scale on the basis of GHG intensi- 
ties at low and deep water-table conditions, which reflect wetland 
drainage/desiccation (WTL-3: <-70 cm; Extended Data Table 1 
and Supplementary Fig. 3). Over the past 71 years (1950-2020), 
46.22% of global wetlands have been degraded (4.85 Mkm/), releas- 
ing 276.4+ 175.5 GtCO,eq (95% confidence interval of GHG emis- 
sions) to the atmosphere (Fig. 3a and Extended Data Table 2). 
Russia, Brazil and Canada were the largest emitters because of 
their vast wetland areas with a high density of soil organic carbon, 
contributing to nearly one-half of global wetland GHG emissions: 
18.6%, 15.1% and 14.6%, respectively (Figs. 3b and 4a and Extended 
Data Fig. 5a). 

Following a history-derived, business-as-usual scenario for the 
future, we project that continued wetland degradation (7.76 Mkm?, 
74.0%) would release an estimated total of 407.9+ 251.5 GtCO,eq 
into the atmosphere from 2021 to the end of the twenty-first cen- 
tury (Fig. 4b). Of these, 155.6 GtCO,eq (38.1%) will be emitted from 
freshwater marsh and floodplain and 96.7 GtCO,eq (23.7%) from 
peatlands, with the latter emitting an average of 1.21 Gtyr™’, which 
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is consistent with 1.32Gtyr! or 1Gtyr' reported by ref. '* and 
ref. **, respectively. Carbon dioxide would contribute the highest 
emissions: 306.1+159.4Gt (Extended Data Fig. 6 and Extended 
Data Table 2). Regionally, 71.1%, 4.6% and 24.3% of the GHG emis- 
sions would be from boreal, temperate and tropical regions, respec- 
tively. The estimate for the tropics is lower than that reported by 
ref. '° because we consider the depletion of the carbon pool during 
2021-2100 (Extended Data Fig. 7). Furthermore, we recognize the 
potential for strong, positive climate feedback in boreal areas stem- 
ming from the loss of substantial carbon storage in the future’. 


Emissions reduction potential under rewetting scenarios 

To explore the potential for reducing GHG release from degraded 
wetlands, we consider two peak clipping schemes based on the 
rewetting of all degraded wetlands (ALL) and rewetting wetlands 
that contain only high organic carbon stocks (high-OCS) (Methods 
and Supplementary Fig. 3). We find that of a total of 4.85 Mkm* of 
degraded wetlands until 2020, fewer than half (2.02 Mkm?) were 
still emitting GHGs, and the remaining 2.83 Mkm? were completely 
degraded (Supplementary Fig. 4). In the case of the latter, the soil 
carbon pool limits the duration for which GHGs are potentially 
emitted from ecosystems***’. We then estimate that a widespread 
rewetting of degraded wetlands, with restoration rates of the same 
magnitude as the historical degradation rates, can potentially 
reduce GHG emissions by 248.7+154.6GtCO,eq (ALL rewet- 
ting scenario) and 156.4+94.2 GtCO,eq (high-OCS scenario) by 
2100 (Fig. 4c and Extended Data Table 2). The latter contribution 
was mainly from freshwater marshes and floodplains (48.6%) and 
peatlands (38.2%). The reduction of emissions from peatlands is 
an average 0.75GtCO,eqyr’, consistent with 0.5Gtyr' reported 
by ref. *! based on the ‘optimal rewetting’ scenario in which 65% 
of peat is under cropland and grassland (0.77Gt in 100% peat). 
The corresponding CO, reductions are 192.9+ 104.4 Gt (account- 
ing for 77.6% of the sum of the three GHGs at 248.7 GtCO,eq) 
and 107.0+48.7Gt (68.4%), respectively, for the two scenarios. 
Although the area proportion of high-OCS to ALL wetlands is only 
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Fig. 3 | GHG emissions from degraded wetlands under different scenarios. 
a, Time series of GHG emissions from degraded wetlands under three 
scenarios since 1950. The emissions are constrained by natural WET 

index and soil organic carbon pool. ‘History trend’ is history-derived 
scenario. ‘Rewet (ALL)' and ‘Rewet (high-OCS)' are based on the rewetting 
restoration of all and only high organic carbon stock wetlands, respectively. 
b, GHG net flux from degraded wetlands in main countries and continents 
over different periods. ‘Others’ refers to the sum of GHGs from countries 
that were not in the top ten of GHG emitters. EU, Europe; NA, North 
America; SA, South America; AS&OA, Asia and Oceania; AF, Africa. 


42.9% (-0.92 +2.92 Mkm’ versus 1.34+2.92 Mkm?, Supplementary 
Fig. 4), the GHG emissions reduction potential is as high as 62.9% 
(156.4 GtCO,eq versus 248.7 GtCO,eq) owing to the differences in 
the higher carbon density. 

These projections are presumed to consider the effects of various 
types of wetlands and a changing climate. The duration of potential 
degradation related to GHG emissions from wetlands is constrained 
by the initial carbon pool and degradation rate (equation (4) and 
Supplementary Fig. 5), which vary according to wetland type and 
the climate regime. However, in our extension of this assessment 
to the future, we assume the effect of climate change will be neg- 
ligible on emission rates, although there would be an unknown 
additional climate effect. We base this assumption on the similar- 
ity in GHG emissions rates between temperate and boreal regions 
for drained conditions (WTL-3: <-70 cm) (Fig. 2a,b and Extended 
Data Table 2). The similarity suggests emissions in cold regions 
would not change greatly with the magnitude of anticipated warm- 
ing. Although the same comparison applied to the temperate versus 


tropical wetlands indicates a non-negligible impact of warming, the 
overall impact should be limited because of the low proportion of 
total emissions from temperate wetlands (4.6%; Extended Data 
Fig. 7). By contrast, emissions of GHGs in the tropics would prob- 
ably be sufficient to deplete the carbon pool before 2100 (Extended 
Data Fig. 7a,d). Another issue of consideration is the effect that 
individual wetland categories would have on the effectiveness of the 
rewetting measures. However, there are consistent GHG exchange 
rates across various wetland categories (Supplementary Fig. 6a-f) 
that are similar to the sum exchange of GHGs across different cli- 
mate regimes (for example, WTLO; Fig. 2a-c and Extended Data 
Table 2), although for which CO, and CH, fluxes vary significantly 
(Supplementary Fig. 7a-c). 

A general limitation in our assessment stems from the natural 
WET index being produced at an intercontinental scale rather than 
at a local scale*’. This limitation prevents the development of a more 
detailed assessment of the GHG budgets associated with wetland 
degradation and/or restoration. In addition, we could not consider 
fire disturbances for wet versus dry conditions’ or the response 
of substratum over one metre deep. Thus, our estimated GHG 
emissions from degraded wetlands may be conservative, and the 
GHG-reduction potential of rewetting programmes is likely to be 
underestimated. Nevertheless, our estimated annual emissions from 
degraded wetlands for CO, and N,O, or from whole natural freshwa- 
ter wetlands for CH,, are significantly correlated with their respec- 
tive annual growth rates in atmospheric concentration” during the 
past three decades (P< 0.05; Extended Data Fig. 8; more details in 
Wetland GHG budgets and inter-annual atmospheric GHG growth 
rates), supporting the utility of our method of integrating the natural 
WET index with the empirical GHG emission rates. Furthermore, 
the magnitudes of the estimated CO,, CH, and N,O emissions from 
the wetlands, equalling 10.8+ 6.2%, 38.5+ 16.7% and 30.5+ 19.4%, 
respectively, of the anthropogenic sources'*"'”’, are in line with the 
aforementioned correlations. 

Despite uncertainties, we find that wetland rewetting is an effec- 
tive nature-based solution to mitigate climate change. The rewet- 
ting ALL scenario and the high-OCS scenario require preserving 
or restoring 4.26 Mkm? and 2 Mkm?’ areas of degraded wetlands, 
respectively, compared with the business-as-usual scenario derived 
from historical trends (Supplementary Fig. 4). The two rewetting 
scenarios can reduce GHG emissions by 583.8tCO,eqha! and 
782 tCO,eq ha“, respectively. The potential of GHG reduction from 
wetland restoration at this scale is higher than that from the reha- 
bilitation of other types of ecosystems, for example, forest regrowth 
equivalent to 394.2tCO,ha™ across 6.78 Mkm? of afforestation 
under the ‘maximum scenario, or 504.3 tCO, ha” across 3.49 Mkm? 
under the ‘national commitments’ scenario**. At present, Indonesia, 
Europe and North America have already shown the benefits of rais- 
ing water tables by both artificial and natural means‘. In cases 
where the water table is lowered by groundwater extraction, water 
conservancy measures may be needed to regulate water use—and 
the potential negative impacts of doing so should be balanced*””. 


Conclusion 

In conclusion, the nonlinear thermal-wetness influence on wet- 
land GHG fluxes, whereby a near-surface water table produces 
near-neutral GHG flux across broad temperature gradients, sug- 
gests that rewetting wetlands is an effective nature-based solution 
to mitigate climate change. A volume equivalent to 10% of anthro- 
pogenic CO, emissions could be reduced through wetland restora- 
tion. By quantifying the impact of natural wetland area changes on 
multiple GHG budgets under several scenarios, we provide primary 
information for nature-based solutions predicated on wetland res- 
toration for countries aiming to achieve net-zero emission targets”. 
Furthermore, we emphasize the enormous loss of organic matter 
and GHG emissions from over half of the global wetland ecosystems 
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Fig. 4 | Spatial pattern of the GHG emissions owing to wetland degradation and reduction potential via rewetting wetlands. a,b, The GHG emissions 
owing to wetland degradation under history-driven scenario in 1950-2020 (a) and under history-driven, business-as-usual scenario in 2021-2100 (b). 
(¢) The reduction potential under Rewet (ALL) scenario in 2021-2100. 
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due to drying. Finally, we show that future emissions from these 
sources can be mitigated or even halved by rewetting wetlands to a 
near-surface water table. 
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Methods 

In situ database of GHGs. We searched for literature through the Web of Science 
using the following strings: (greenhouse gases flux (greenhouse gas* OR GHG* 
OR carbon dioxide OR CO, OR net ecosystem productivity* OR NEP OR net 
ecosystem exchange* OR NEE OR carbon OR methane OR CH, OR nitrous 
oxide OR N,O OR flux* OR emission* OR global warming potential OR GWP) 
OR (degradation* OR decline*) and wetland (wetland* OR water* OR peatland* 
OR bog* OR fen* OR swamp* OR mire* OR soil* OR river* OR paddy OR 

pool* OR floodplain* OR reservoir* OR coastal* OR saltmarsh*)). The search 
returned 33,835 papers after excluding those in irrelevant fields. After screening 
the manuscripts to ensure that the records contained measured data from field 
monitoring projects, 2,563 papers were manually selected, from which we read to 
extract the following information: (1) gas sample collection methods (this content 
may be missing, for example, eddy covariance method (EC) measured CO, and/or 
CH, fluxes without collecting gas samples) and measurement techniques; (2) time 
intervals and the duration of sample collection (we specified that the monitoring 
span during the growing season should be at least three months in temperate 
regions); (3) detailed site information, in particular geographic coordinates or 
maps that could be georeferenced, as well as topography, ecosystem type, plant 
species and soil carbon features. After screening, we obtained 504 papers with valid 
data (Supplementary Data 1). Parts of those data were extracted from figures by 
GetData Graph (version 2.26). 


Record details. We constructed a GHG net-flux database consisting of 3,704 
site-years (1,875 sites) for locations situated between 81°48’N and 78°01’S 
(Extended Data Fig. la—c) across all the continents, but mainly in Asia, Europe 
and North America. The sample collection method and measurement technique 
used in these field monitoring projects usually depended on the type of wetland 
and associated environmental conditions. Common collection methods included 
the static chamber method”', dynamic chamber method” and floating chamber 
method*’. Common measurement techniques involve the use of infrared gas 
analyser and the combination of gas chromatography with physical model methods 
(for example, EC™, CO,/CH,-diffusivity formula™, oxygen diffusivity formula” or 
chlorophyll-dissolved oxygen model”). We assume that these GHG exchange rates 
from peer-reviewed papers based on different sampling methods are of equivalent 
accuracy, although there may be differences in precision. 

Nomenclature corresponding to the exchange of CO, between the ecosystem 
and the atmosphere used in various methods is diverse. The EC method generally 
uses a net ecosystem exchange (NEE) to characterize the flux, while the chamber 
method adopts the net ecosystem production approach. While a few sites measure 
both vertical exchange and transverse flow of GHGs**’, most of the sites report 
only the measurements of vertical fluxes. Consequently, the datasets we built 
describe the fluxes in the vertical direction. In cases where data were reported for 
the growing seasons (Extended Data Fig. 1b), we filled in the missing data using 
linear regression on the basis of complete data containing both growing and annual 
records (Supplementary Table 1). There was no significant interaction effect of 
climate regimes on the slope, and only a subtle effect on the intercept of the NEE 
regression equation for temperate zones. Therefore, we chose linear regression 
equations that did not consider the interaction of the climate regimes. 

We used NEE to describe the net vertical exchange capacity: 


NEE = Re — GPP = Rh — NPP (1) 


where Re is ecosystem respiration, GPP is gross primary productivity, Rh is 
heterotrophic respiration and NPP is net primary productivity. Positive values of 
NEE indicate carbon loss via CO, emission to the atmosphere, whereas negative 
values indicate carbon gain because CO, is retained/stored in the ecosystem. 

Building on equation (1), we calculated the total flux of three GHGs—carbon 
dioxide, methane and nitrous oxide—as 


GHGs = NEE + 34CH, + 298N20 (2) 


where CO, is substituted by NEE. Units for all terms are kgCO,eq ha" yr’; the 

weights 34 and 298 are global warming potential (GWP) for CH, and N,O to 

CO, eq by weights on a 100 yr perspective with feedbacks considered”, respectively. 
Uncertainty is represented by incorporating 95% confidence intervals (CIs): 


GHGs_CI = NEE_CI + 34CH4_CI + 298N20_CI (3) 


where the Cl of three GHGs on the right side of equation (3) originate from within 
groups, whereby uncertainties in empirical GHG emission rates are obtained for 18 
groups, including 6 WTL categories x 3 climatic zones. 


Wetland categories. The majority of the records collected in this study were 
measured in natural wetlands, including bogs, fens, mires, swamps, marshes and 
floodplains; we also considered water bodies as extensions. We distinguished 
peatlands from non-peatlands on the basis of the soil organic layer thickness 
(peatlands: >40 cm)” to explore the influence of soil organic matter on multi-GHG 
fluxes. We divided all sites into six water categories on the basis of either the 

water level (flooded) or the water-table height (non-flooded) relative to the 
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surface (WTL) for three climatic regions (tropical, temperate, and boreal or 
high-altitude areas, using 4°C and 17°C as thresholds of multiyear-average annual 
temperature)”. We also treat water level as a continuous variable in Extended Data 
Fig. 2, using 2,318 site-year records that reported exact water levels. Unfortunately, 
reports with high-quality water-level records from the tropics accounted for only 
7% of the total. As such, the trend for the tropics using the continuous variable 
approach has more uncertainty than for other climates. Therefore, we segmented 
all the WTL data. The WTL classes were separated by mean water table/level 
during the growing season due to the discontinuity of the reported hydrological 
records during winter at most sites. Six classes are designated: WTL-3, <-70 cm; 
WTL-2, -70cm < WTL-2 <-50cm; WTL-1, -50cm < WTL-1 <-30cm; 

WTLO, -30cm <WTLO <-5 cm; WTLI, -5cm<WTLI1 < 40cm; and WTL2, 
>40cm. Positive numbers indicate that the water level is above the surface. This 
classification is based on the empirical water-level critical value of the soil moisture 
conditions, ranging from drought to moist to near-saturated to oversaturated®. We 
used an empirical value of 40cm water level as a threshold to classify high water 
levels as they usually inhibit the establishment of emergent plants. We defined the 
start/end of a growing season as the time when the daily mean temperature for five 
continuous days was above/below 5°C for the first time". 

We extracted the global natural wetland map from the GLWD level 3 
(GLWD-3, classes 4, 5, 8-12)*, which provides reliable areas of global wetlands, 
with tropical peatlands distinguished. Its classification of wetlands also agrees 
with most records in the literature we selected. Gridded long-term mean annual 
air temperature was calculated using the monthly data from the latest European 
Centre for Medium-Range Weather Forecasts reanalysis (ERA5 (ref. °°)) for 1978- 
2018. Gridded wetland soil OCS was obtained from SoilGrids”’, which indicates 
that natural freshwater wetlands cover 7.66 million km? and store 329.5 GtC in the 
upper one metre of the soil (Extended Data Table 2). The proportion of wetlands in 
boreal, temperate and tropical regions are 45.8%, 10.8% and 43.4%, respectively. 


Dynamic area of wetlands under three scenarios. Means and 95% CIs of GHG 
net fluxes in WTL-3 (<-70 cm) for different climatic zones were used as the 
emission potential from degraded wetlands (Extended Data Fig. 2b and Extended 
Data Table 1). We employed the mean OCS for each wetland category in each 
country as a constraint and then calculated the duration potential (Supplementary 
Fig. 5) of soil carbon efflux from degraded wetlands for 21 wetland types (3 climate 
zones (tropical, temperate, and boreal or high-altitude zones) x 7 wetland category 
groups (bog, fen, mire, swamp forest, flooded forest, freshwater marsh and 
floodplain)) at the national scale as follows: 


DP = (OCS — OCS,) /NEE, DP > 1 (4) 


We used an OCS, (threshold) of 50 tha” (ref. °°), after which an ecosystem has 
no potential for net CO, emission from the soil layer. To avoid double counting, 
we did not consider the other GHG exchanges of completely degraded wetlands, 
which have been definitely converted into farmland or pasture, and so on. These 
are included in the estimate of land-use and agriculture emissions by the Carbon 
Budget Project** and the Food and Agriculture Organization”. 

We rebuilt the historical trend and predicted the future wetland degradation 
rate at the continental scale using the natural WET index" for 1970-2015 
(Supplementary Fig. 2). The WET index is a multi-source composite index 
to represent the proportion of wetland degradation or construction during 
1970-2015 (Europe in 1970-2013). We extended the time series to 1950-1969 and 
2016-2100 from regressions using the 1970-1990 and 2000-2015 data (Europe 
predicted by 2000-2013 due to the data restriction); these calculations are for 
the history-derived (history-derived) scenario. We used two wetland restoration 
benchmarks: one considered all wetlands (ALL) and the other involved only the 
high-OCS wetlands (high-OCS) (Supplementary Fig. 3). High-OCS wetlands were 
determined as those with a duration potential >80 yr. The total area of high-OCS 
wetlands worldwide is 3.29 Mkm’, of which 34.5% are degraded. The ALL and 
high-OCS scenarios were both grounded on the assumption that the restoration 
rate of those degraded wetlands would exceed that of the historical degradation 
rate in all continents by 2030. Thus, before 2030, restoration proceeded at a rate 
equivalent to the absolute value of the degradation rate during 2000-2015 (2000- 
2013 for Europe). The increase stopped when the natural WET index recovered to 
the level in 1950 (Supplementary Fig. 2). 

A completely degraded wetland was defined as having no potential for soil 
carbon loss, and the wetland during degradation was defined as having continuous 
loss of carbon and nitrogen (WTL-3) to the atmosphere. The remainder were 
classified as initial and rewetted wetlands. The component of dynamic wetland area 
under three scenarios in Supplementary Fig. 4 shows the trends of the area in three 
conditions during 1950-2100 (more details in Supplementary Data 1). 


Supporting calculations. Empirical GHG exchange rates. On the basis of the 
empirical values (WTLO, WTL2) in the relationship between WTL and GHG 
emissions across different temperature regimes (Extended Data Table 1), we 
calculated that the GHG emissions from global water bodies (lakes and reservoirs) 
are 1.0+0.93 GtCO, and 127.5+49.4TgCH,, and CH, emissions from natural 
freshwater wetlands are 144.4+67.9Tg. These results agree with previous reports, 
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whose corresponding values are 1.2*};3! GtCO, (ref. **), 175.2481 TgCH, (ref. **) 
and 148.6 + 15.2 TgCH, (ref. '*), respectively (Supplementary Fig. 1). 

The GHG emissions from boreal and temperate regimes in WTL-3 
(<-70cm) are 19.7 and 11.2tCO,eqha™ yr", respectively, which are similar to 
the drained-induced GHG emission factor 16.1 tCO,eq ha" yr" for boreal and 
temperate in ref. '’. Furthermore, the CO, emission rate in the boreal regime is 
1 3.43128 tCO, ha" yr“, which is similar to the cultivated northern peatlands 
emission factor 13.2+°73tCO,ha™ yr” in ref. ®’ from a process-based land surface 
model. 


Estimation of GHG emissions from degraded wetlands under three scenarios. We 
estimated the changes of GHG emissions driven by the degraded wetland area and 
OCS in various wetland types under three scenarios (history-derived, ALL and 
high-OCS) and across different scales. At the national scale, countries with large 
OCS were the dominant emitters of GHGs, mainly CO,. The top ten GHG emitting 
countries contributed to 79.6% of the total emissions in the period 1950-2020 
(Fig. 3 and Extended Data Figs. 5 and 6). At the continental scale, the top two 
continents, South America (SA) and North America (NA), accounted for 37.7% 
and 22.7%, respectively. With respect to climate regime, countries in tropical and 
boreal regions occupied 55.4% and 39.7%, respectively, due to their potent outflow 
and high OCS. 

In the scenario for which the historical trend continues during the period 
2021-2100 (history-derived scenario), the wetland degradation area will increase 
to 74.0% by the end of the twenty-first century, and the induced GHGs will be 
enlarged by 1.48 times to 407.9 Gt (Figs. 3 and 4 and Extended Data Table 2). 

In addition, the present geographical pattern of outflow will change: in tropical 
regions, along with the loss of most soil OCS, emissions will reduce from 155.2 Gt 
to 99.1 Gt, while boreal regions will become the dominant emissions source, 
increasing from 107.9 Gt to 290.0 Gt (39.0% to 71.1%). 

Under the wetland restoration scenarios for which rewetting occurs for all 
degraded wetlands (ALL) or only high-OCS degraded wetlands (high-OCS), the 
total GHG emissions could reduce by 248.7 or 156.4 Gt, respectively. The emissions 
reduction results mainly from the reduced GHG emissions over boreal regions. 

In the high-OCS scenario, CO, accounted for 79.2% of the global total GHG 
emissions (Extended Data Table 2), higher than in the ALL scenario (71.2%) and 
the history-derived scenario (75.0%). This finding results from a large percentage 
of the degraded low-OCS wetlands being distributed in the tropics but not restored 
in the high-OCS scenario, where the tropics contribute to the largest portion of 
global CO, emissions with the highest CO, outflow potential per area (Fig. 2 and 
Extended Data Fig. 5). 

All the preceding conclusions are based on a GWP for CH, of 34-fold that of 
CO, by weights over a 100 yr period. We have also supplemented our assessment 
of emissions and reductions with a 28-fold GWP or a 45-fold sustained GWP*. 
The results of GHG emissions estimation show differences of +0.23% and —0.43%, 
respectively, compared with the history-derived scenario with a 34-fold base in 
1950-2100. The estimation of GHG emissions reduction change is +0.37% and 
-0.67% under the rewetting ALL scenario, and +0.22% and -0.41% under the 
rewetting high-OCS scenario in 2020-2100, respectively. 


Wetland GHG budgets and inter-annual atmospheric GHG growth rates. There 

are significant correlations between annual growth rates in atmospheric 
concentrations and changes in the wetland-induced flux for CO,, CH, and N,O in 
1979-2018 (P<0.05; Extended Data Fig. 8a-i; N,O is in 1979-2016). The sum of 
N,O emissions from degraded wetland and Food and Agriculture Organization 
agriculture-total” is also significantly correlated with atmospheric concentration 
growth rates (P< 0.05; note that N,O data exclude two early extreme values in 1979 
and 1982). Emissions of CO, from degraded wetlands are highly consistent with 
those from land-use changes in the Carbon Budget Project’, with the former being 
~81.4% (ranging from 64.5% to 97.5%) of the latter across 30 years. The net CH, 
emissions from degraded and from initial and rewetted wetlands both exhibited 

a downward trend during 1979-2018. This change may have contributed to the 
decline of atmospheric methane growth rates before 2005, which is reversed by the 
increase of emissions from other major sources (for example, agriculture”) since 
then. The evident correlations between wetland budgets and atmospheric growth 
of three key GHGs indicate the non-negligible impact of wetland degradation. 
Indeed, CO,, CH, and N,O emissions from wetlands were equal to 10.8 + 6.2%, 
38.5 + 16.7% and 30.5 + 19.4% of those from anthropogenic sources'*“!’, similar 
to their contributions of 19%, 21% and 34% of the variation in atmospheric 
concentrations during 1979-2018 (Extended Data Fig. 8). 


Uncertainties. Wetland area. Although many global wetland area products using 
diverse classification rules have been released**”””*, uncertainty remains in wetland 
characterizations and distribution worldwide. In particular, human-made wetlands 
have been increasing greatly in recent years, and natural wetlands continue to 
degrade**”°. The natural wetland classes (bog, fen, mire, swamp/flooded forest, 
freshwater marsh and floodplain) of GLWD-3 were built from the following 

three datasets in the 1990s: ArcWorld”, Digital Chart of the World’* and World 
Concervation Monitoring Centre”. We conservatively considered the area of 
wetlands in GLWD-3 as the background value in 1990. The GLWD-3 did not 


include wetlands smaller than 0.1 km?; however, this threshold ensured that the 
hydrologic features and biogeochemical processes of degraded wetlands can be 
restored to their initial states in a short term*’. 


Degradation trend of natural wetlands. Many countries lack baseline wetland 
inventories that allow us to accurately track the lengthy and complex degradation 
of natural wetlands". The WET index represents area change based on over 2,000 
wetland area records in long-term time series from six regions and distinguishes 
between human-made and natural wetland changes from 1970 to 2015. Given 
the differences between natural and artificial wetlands in the basal features (for 
example, OCS and hydrologic features) and disturbance (for example, artificial 
landscape and artificial nitrogen input), we considered only natural wetlands via 
down-scaling the natural WET index from six continents to countries and/or 
regions. Therefore, additional surveys and remote-sensing data would provide a 
more accurate assessment at a finer scale in the future. 

According to the WET index, the degraded wetland areas were 4.85 Mkm? 
for 1950-2020 (46.22% of the global natural wetlands). Under a history-derived, 
business-as-usual scenario for the future, we projected that continued wetland 
degradation will reach 7.76 Mkm? (74.0%) by the end of the 2021-2100 period. 
Note that these overall estimates include those wetlands degraded not only by 
land-use change (for example, reclamation and draining) but also by other factors 
affecting the hydrological characteristics of wetlands (desiccation). For example, 
the peat loss to extraction or farming in Europe accounts for only 11% of European 
peatlands”, but the water levels in another 50% of European peatlands are also 
declining, causing degradation*’. 


Wetland categories. Empirical parameter generalization is based on the premise 
that there is no difference in emission potential across various wetland categories. 
Indeed, we found that, in the WTLO group, where GHG emissions are close to 
neutral (Fig. 2 and Extended Data Table 1), there is almost no significant difference 
in GHG emissions among the main wetland categories for both tropical and 
temperate regimes (Supplementary Fig. 6). However, because of the lack of various 
types of wetlands in boreal regions (dominated by peatlands), we did not test the 
differences in GHG emissions across diverse categories for boreal climates. In the 
WTL-3 group, due to the lack of various wetland categories for comparison, we 
used the OCS to calculate the duration potential, which is further used to constrain 
the empirical parameter generalization for the estimate of the potential GHG 
emissions from degraded wetlands (equation (4)). 


Data availability 

GLDW dataset is available at http://www.wwfus.org/science/data.cfm. 

Soilgrids dataset is available at https://soilgrids.org. ECMWF reanalysis climate 
data are available at https://cds.climate.copernicus.eu/#!/home. FAOSTAT 
emissions database is available at http://www.fao.org/faostat/en/#data/GT. 
Atmospheric concentrations data are available at https://ourworldindata.org/ 
atmospheric-concentrations. All GHG data are available in the main text or 

the supplementary materials. The database of global, in situ, GHG exchange 
information for wetlands, drawn from 3,704 site-year records, is summarized in 
Supplementary Data 1. Source data are provided with this paper. 


Code availability 

The scripts used to generate all the results are MATLAB (R2018a), R-4.1.0 and 
Python 2.7 based on arcpy. Analysis scripts are available at https://github.com/ 
XiaoBai0417/Multi-greenhouse-gas-assessments. 
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Extended Data Fig. 2 | Annual GHG flux values versus water level in wetlands. a, GHG balance, the nonlinear trend is formed by 103 site-year records 
reporting exact water levels and complete data including three major GHGs. b, CO,, CH, and N,O exchange. NEE (e), CH, (d) and N,O (e) exchange across 
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Extended Data Fig. 5 | Greenhouse gas emissions from degraded wetlands in countries (a) and wetland categories (b). a, The country’s historical 
emissions. The color of each circle corresponds to the axis of the same color (red/right; blue/left). The size of a circle represents the amount of soil 
organic carbon stock. b, Emissions from different wetland types under three scenarios. For details see Supplementary Data 1. 
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Extended Data Fig. 8 | Inter-annual atmospheric GHG concentration changes and emissions from wetlands and other major sources. a-i Data for CO, 
(a-c), CH, (d-f) and N,O (g-i). Significant correlations exist between the flux of CO,, CH, and N,O from wetlands and their respective atmospheric 
concentration growth rates in the past three decades. In addition, the correlations persist when adding other major emission sources to wetlands. Note 
that atmospheric growth rates of N,O for 1979 and 1982 are excluded. 
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Extended Data Table 1 | Wetland greenhouse gas (GHG) net fluxes in different climate regimes under various water table levels 
(WTL) 


Boreal Temperate Tropical 


WTL GHGs 95%CI GHGs 95% CI GHGs 95%CI 


2 8.89 4.82 2o.02 10.45 34.90 18.78 
1 6.41 2.24 9.39 5.56 8.89 16.89 
0 5.10 4.36 4.27 2.46 10.90 14.37 
-1 7.45 3.48 12.21 4.69 25:20 8.20 
-2 11.69 9.38 i937 11.45 13.07 19.39 
-3 19.67 11.34 11.25 13.18 61.20" 38.50" 


unit: tCO,eq ha+ yr+; * data from J. Leifeld et al., 2019 (ref. 19). 
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Extended Data Table 2 | Wetland characteristics and GHG emissions for each continent 


Continent 


Europe 
North America 


Latin America 


Asia & 
Oceania 


Africa 
Total 


Reduction 


Area 
Gnillion km:) 
1.27 
2.71 
1.37 
1.23 
1.08 


7.66 


ocs* 
(GtC) 
99.51 
158.81 
31.05 
26.88 
13.23 


329.48 


1950-2020 
discharged 


53.94 + 32.66 
(68.38%)t 
62.64 + 40.67 
(70.19%) 
104.15 + 65.83 
(95.62%) 
36.07 + 23.96 
(89.83%) 
19.64 + 12.37 
(95.73%) 
276.44 + 175.49 
(83.80%) 
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GHG net-flux from degraded wetland (Gt CO.eq) 


2021-2100 
projected from 
historical trend 
199.70 + 120.69 

(68.38%) 
100.77 + 61.55 

(68.69%) 
52.39 + 33.17 

(95.57%) 
35.70 + 23.93 

(89.03%) 

19.30 + 12.18 

(95.67%) 
407.87 + 251.51 

(75.05%) 


2021-2100 2021-2100 
rewetting all wetlands _ rewetting high OCS 
scenario wetlands scenario 
66.11 + 39.93 66.40 + 40.26 
(68.37%) (68.39%) 
75.05 + 45.28 78.25 + 48.36 
(68.41%) (68.79%) 
9.95 +6.31 52.39 + 33.17 
(95.50%) (95.57%) 
6.17 + 4.22 35.08 + 23.32 
(87.34%) (89.38%) 
1.87+ 1.18 19.30 + 12.18 
(95.52%) (95.67%) 
159.14 + 96.92 251.42 + 157.28 
(71.14%) (79.20%) 
248.72 + 154.59 156.44 + 94.23 
(77.55%) (68.37%) 


*: OCS represents organic carbon stocks in soil layers from zero to one meter deep. *: The figures in the parentheses indicate the percentage of CO, 


emissions. 


The figures in the parentheses indicate the percentage of CO, emissions, OCS represents organic carbon stocks in soil layers from zero to one meter deep. 
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